% 给定参数
V = 10; % 船舶速度（米/秒）
t = 10; % 转弯/跳投时间（秒）
theta = pi/6; % 转弯角/跳投角（弧度）
m = 10000; % 船舶排水量（千克）
C_D = 0.5; % 阻力系数（无量纲）
A = 10; % 船舶的迎风面积（平方米）

% 计算转弯半径
R = V * t / (9.81 * tan(theta));

% 计算跳投轨迹
time = linspace(0, t, 100); % 时间点
x = V * time * cos(theta) - 0.5 * C_D * A * V^2 * time.^2 / (m * 9.81);
y = V * time * sin(theta) - 0.5 * 9.81 * time.^2;

% 绘制轨迹图
figure;
subplot(2, 1, 1);
plot(time, x, 'b-');
xlabel('Time (seconds)');
ylabel('Horizontal Position (meters)');
title('Horizontal Trajectory');

subplot(2, 1, 2);
plot(time, y, 'b-');
xlabel('Time (seconds)');
ylabel('Vertical Position (meters)');
title('Vertical Trajectory');